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Gusts  from  boundary-layer-scale  convective  eddies  modify  the  structure 
of  the  atmospheric  surface  layer  by  intermittently  intensifying  or  dimin¬ 
ishing  the  local  wind  speed.  A  simple  model  for  the  effects  of  these  gusts 
is  proposed,  based  on  the  following  two  assumptions:  (1)  the  wind  gusts 
have  an  isotropic  Gaussian  probability  distribution  with  standard  devia¬ 
tion  proportional  to  re*  (the  convective  velocity  scale)  and  (2)  the  surface- 
layer  wind  and  temperature  profiles  attain  local  equilibrium  with  the  wind 
gusts.  The  minimum  friction  velocity  predicted  by  the  model  has  the  same 
dependence  on  surface  roughness  predicted  by  Schumann's  earlier  "slab" 
model  for  convective  boundary  layers.  However,  the  current  model  also 
applies  to  situations  where  the  mean  wind  is  nonzero.  It  predicts  the  break¬ 
down  of  global  Monin-Obukhov  similarity  for  the  surface-layer  wind  shear 
and  temperature  gradient  in  highly  convective  conditions  {Ur/w^  approxi¬ 
mately  1  or  smaller,  where  Ur  is  the  mean  wind  speed  at  the  top  of  the 
surface  layer).  Also  in  contrast  to  existing  similarity  theories,  the  horizon¬ 
tal  wind  variance  exhibits  a  significant  dependence  on  height  and  surface 
roughness  near  the  ground,  even  for  moderate  convection.  The  tempera¬ 
ture  variance  is  nearly  unaffected  by  the  gusts,  because  of  its  weak  depen¬ 
dence  on  the  local  wind  speed  in  convective  conditions. 
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1.  Introduction 


One  of  the  primary  goals  of  micrometeorology  is  to  develop  statistical  de¬ 
scriptions  for  the  surface-layer  (near-ground)  structure  of  the  atmosphere. 
These  statistical  descriptions  have  many  valuable  applications,  such  as  fore¬ 
casting  the  turbulent  transport  of  biological  and  chemical  agents  and  cal¬ 
culating  the  propagation  of  acoustic  and  electromagnetic  waves.  The  most 
successful  model  for  surface-layer  statistics  is  the  turbulence  similarity  the¬ 
ory  developed  by  Monin  and  Obukhov  (1954),  which  hypothesizes  that  the 
statistics  depend  on  just  four  parameters:  the  friction  velocity  u*,  the  sur¬ 
face  temperature  flux  Qs,  the  buoyancy  coefficient  b  =  g/Ts  (where  g  is 
gravitational  acceleration  and  Tg,  the  surface  temperature),  and  the  height 
from  the  ground  z.  Note  that  no  parameters  in  the  Monin-Obukhov  sim¬ 
ilarity  theory  are  related  to  the  large-scale  structure  of  the  atmospheric 
boundary  layer.*  The  theory  assumes  that  surface-layer  processes  are  not 
significantly  affected  by  the  large  boundary-layer-scale  eddies  (also  called 
thermals  or  -plumes)  characteristic  of  convective  conditions.! 

Primarily  during  the  1960s  and  1970s,  semiempirical  equations  for  the  ver¬ 
tical  profiles  of  mean  wind  and  temperature  were  developed  on  the  basis 
of  Monin-Obukhov  similarity.  These  equations  fit  the  vertical  profiles  quite 
well  and  have  since  come  into  widespread  use.  Other  statistics,  such  as  the 
temperature  and  vertical  velocity,  have  also  been  found  to  agree  with  the 
theory.  Nonetheless,  experiments  have  clearly  demonstrated  that  not  all 
near-ground  statistics  obey  Monin-Obukhov  similarity.  Most  significantly, 
the  theory  does  not  work  well  for  the  variance  and  spectrum  of  the  horizon¬ 
tal  wind  components,  which  have  been  found  to  depend  on  Zi,  the  height 
of  the  boundary-layer  capping  temperature  inversion  (Kaimal  et  al,  1976; 
Caughey  and  Palmer,  1979;  Khanna  and  Brasseur,  1997).  Furthermore,  Bel- 
jaars  (1994)  and  Grachev  et  al  (1998)  demonstrated  that  parameterizations 
for  surface  fluxes  and  drag  coefficients  in  convective  conditions  agree  bet¬ 
ter  with  observations  when  the  wind  speed  is  replaced  by  a  scalar  average 
proportional  to  the  mixed-layer  velocity  scale,  tc*  =  (QsZig/Tg)^^^,  which 
is  representative  of  the  large  eddies. 

Potential  limitations  of  Monin-Obukhov  similarity  in  convective  condi¬ 
tions  were  actually  anticipated  decades  ago,  apparently  first  by  Businger 

*The  boundary  layer  is  defined  as  "the  layer  of  air  directly  above  the  Earth's  surface  in 
which  the  effects  of  the  surface  (friction,  heating,  and  cooling)  are  felt  directly  on  time  scales 
less  than  a  day"  (Garratt,  1992,  p.  1).  It  is  typically  500  to  3000  m  deep  during  the  daytime. 
The  surface  layer  can  be  thought  of  as  the  lowermost  tenth  of  the  atmospheric  boundary 
layer. 

^Convective  conditions  here  means  that  boundary-layer  turbulence  is  being  produced  pri¬ 
marily  by  buoyant  instabilities,  as  opposed  to  wind  shear.  The  buoyant  instabilities  usually 
dominate  when  there  is  solar  heating  of  the  ground. 
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(1973a,b).  The  main  purposes  of  this  report  are  to  systematically  address 
the  behavior  of  surface-layer  statistics  in  convective  conditions  and  to  pro¬ 
vide  a  quantitative  assessment  of  the  breakdown  of  Monin-Obukhov  sim¬ 
ilarity  in  these  conditions.  The  framework  I  adopt  is  an  extension  of  con¬ 
cepts  developed  by  several  previous  authors,  including  Schumarm  (1988), 
Sykes  et  al  (1993),  Beljaars  (1994),  and  Grachev  et  al  (1997,  1998).  The  un¬ 
derlying  principle  is  that  large-scale  convective  eddies  (circulations),  which 
span  the  depth  of  the  boundary  layer,  affect  near-surface  transfer  processes 
by  inducing  a  fluctuating  local  surface  stress.  Grachev  et  al  (1997)  aptly  call 
this  the  convection-induced  stress  regime  to  differentiate  it  from  the  situation 
where  stress  is  induced  primarily  by  a  steady  wind.  To  quote  those  authors: 

Glose  to  the  surface,  a  large-scale  convective  eddy  has  an  in¬ 
stantaneous  local  wind  profile  which  rapidly  achieves  a  steady 
state.  The  local  structure  of  such  a  profile  should  be  essentially 
the  same  as  in  the  usual  mean  wind  profile  (Kraichnan,  1962; 
Businger,  1973a).  As  a  result,  free  convection  can  be  treated  as  a 
particular  case  of  forced  convection  (Sykes  et  al,  1993;  Beljaars, 

1994). 

Therefore  it  is  reasonable  to  apply  the  Monin-Obukhov  similarity  theory 
locally  to  statistics  in  the  near-surface  region  of  the  large-scale  convective 
eddies,  even  though  application  to  global  (ensemble)  average  statistics  may 
be  tenuous. 

In  section  2, 1  show  how  the  premise  of  local  equilibrium,  when  combined 
with  a  probability  distribution  for  the  wind  gusts  induced  by  the  large- 
scale  eddies,  leads  to  a  simple  and  workable  methodology  for  calculating 
surface-layer  statistics.  This  model  is  applied  to  the  calculation  of  mean 
vertical  profiles  and  variances  in  section  3.  Two  appendices  follow  the  body 
of  the  report:  appendix  A  discusses  the  value  of  numerical  constants  used 
in  the  model,  and  appendix  B  provides  approximate  analytical  results  for 
large-scale  eddy  effects  on  atmospheric  statistics. 
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2.  Modeling  Surface-Layer  Gusts 


2.1  Concept 

In  the  conventional  decomposition  of  turbulence  first  proposed  by 
Reynolds,  fields  such  as  wind  and  temperature  are  written  as  the  sum  of 
mean  and  fluctuating  parts.  For  the  wind  velocity  component  aligned  with 
the  mean  flow  direction,  the  Reynolds  decomposition  is 

U  =  {U)+u,  (1) 

where  U  is  the  instantaneous  wind  speed,  u  the  fluctuation,  and  angle 
brackets  indicate  the  ensemble  mean.  If  we  hypothesize  that  the  turbulence 
attains  a  local  equilibrium  with  the  large-scale  convective  eddies,  an  appro¬ 
priate  alternative  decomposition  would  be 

U  =  {U)l  +  u,  (2) 

where  the  subscript  L  on  the  angle  brackets  indicates  a  local  mean  value. 
The  meaning  of  local  here  is  admittedly  imprecise.  The  general  idea  is  that 
a  spatial  filter  is  applied  in  the  horizontal  plane,  eliminating  motions  with 
scales  smaller  than  those  of  the  large-scale  convective  eddies.  A  filter  hav¬ 
ing  a  cutoff  in  the  range  O.lzj  to  0.3zi  would  serve  the  desired  purpose.*  A 
filter  in  time  could  serve  in  place  of  the  spatial  one. 

Since  (u)  is  zero  by  hypothesis,  if  we  take  the  global  average  (ensemble 
mean)  of  equation  (2),  it  follows  that  the  global  average  of  locally  averaged 
quantities  is  simply  the  global  average: 

(u)  =  m^)  ■  (3) 

Throughout  the  rest  of  this  report  I  refer  to  the  variations  in  wind  velocity 
resulting  from  large-scale  convective  eddies  simply  as  gusts.  This  terminol¬ 
ogy  is  borrowed  from  previous  authors,  including  Godfrey  and  Beljaars 
(1991)  and  Grachev  et  al  (1998).  Golloquially,  "gusts"  can  refer  to  variations 
in  the  wind  velocity  other  than  those  originating  from  large-scale  convec¬ 
tive  eddies,  such  as  activity  associated  with  a  thunderstorm  or  a  frontal 
passage. 

Figure  1  illustrates  the  concept  that  local  vertical  profiles  of  wind  and  tem¬ 
perature  obtain  equilibrium  with  gusts,  as  opposed  to  the  actual  mean 

*The  concept  of  applying  a  filter  to  separate  large-scale  motions  from  smaller  ones  is 
discussed  in  more  detail  by  Tong  et  al  (1998).  The  purpose  of  the  filter  in  this  report  is  to 
separate  the  large-scale  convective  eddies  from  other  smaller  scale  motions.  The  purpose 
of  Tong  et  al  was  somewhat  different:  to  mimic  the  resolvable- /subgrid-scale  partition  in 
large-eddy  simulation. 


3 


(1) 


(2) 


(3) 


(4) 


Figure  1.  Vertical 
profiles  of  wind  and 
femperafure  measured 
af  differerenf 
h5^ofhefical  sfafions 
(1-4),  separafed  in 
eifher  space  or  fime. 
Profiles  are  averaged 
with  local  filter.  Upper: 
Profiles  for  a  flow 
having  no  significanf 
gusfiness.  Lower: 
Profiles  for  a  flow 
having  pronounced 
gusfiness.  Local 
femperafure  profiles  are 
in  equilibrium  with 
gusts.  (Arrows 
represent  horizontal 
wind  velocity.  Gray 
scale  represents 
temperature,  with  light 
shading  being  warmer 
air.) 


(1) 


(2) 


(3) 


(4) 


wind,  ([/).  Shown  in  the  figure  are  the  vertical  profiles  of  wind  and  tem¬ 
perature,  after  applying  the  local  filter  {)^,  at  different  measurement  lo¬ 
cations.  For  the  hypothetical  situation  illustrated  in  the  upper  part  of  the 
figure,  the  flow  structure  at  scales  greater  than  the  filter  cutoff  is  negligi¬ 
ble,  so  all  the  observed  local  profiles  are  nearly  identical.  The  wind  speed 
increases  roughly  logarithmically  with  height,  and  a  layer  of  warm  air  lies 
next  to  the  ground.  The  local  profiles  at  each  station  are  in  equilibrium  with 
the  mean  wind  speed.  For  the  situation  illustrated  in  the  lower  part  of  the 
figure,  gusting  is  significant,  and  the  local  wind  and  temperature  profiles 
come  into  equilibrium  with  the  gusts.  At  station  (1),  a  locally  strong  gust 
is  mixing  warm  air  upward  from  the  ground.  At  station  (2),  little  deviation 
of  (C/)^  from  {U)  occurs  so  that  the  wind  and  temperature  profiles  are  the 
same  as  in  the  upper  figure.  At  station  (3),  motion  from  a  large-scale  eddy 
is  in  the  opposite  direction  of  the  prevailing  wind,  hence  decreasing  the  lo¬ 
cal  wind  speed.  As  a  result,  turbulent  mixing  is  weakened,  and  most  of  the 
warm  air  remains  near  the  ground.  If  the  motion  from  a  large-scale  eddy 
is  strong  enough,  the  wind  may  actually  shift  direction  from  the  prevailing 
one,  as  shown  for  station  (4). 

Strictly  speaking,  neither  of  the  situations  illustrated  in  figure  1  pertains  to 
the  customary  application  of  Monin-Obukhov  similarity  theory;  that  the¬ 
ory  applies  to  the  ensemble  flow  characteristics.  Even  though  the  wind  may 
vary  locally  because  of  individual  gusts,  the  flow  is  not  required  to  obtain  a 
local  equilibrium  with  them.  Local  states  of  the  flow  are  just  samples  of  the 
global  ensemble,  and  the  similarity  theory  is  intended  only  for  the  global 
ensemble.  However,  one  might  anticipate  that  Monin-Obukhov  similarity 
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cannot  be  applied  satisfactorily  to  a  global  ensemble  unless  the  local  flow 
samples  are  all  nearly  the  same  as  the  global  average,  as  in  the  upper  part 
of  figure  1.  Otherwise,  the  premise  that  large-  {zr)  scale  eddies  do  not  sig¬ 
nificantly  affect  the  near-ground  flow  is  violated. 

What  justification  does  one  have,  though,  for  applying  Monin-Obukhov 
similarity  to  locally  averaged  flow  fields?  A  simple  supportive  argument 
follows  from  the  relative  time  scales  of  the  large-scale  convective  eddies 
and  the  smaller  surface-layer  ones.  The  convective  eddies  have  a  time  scale 
of  Zijw^.  Taking  Zi  ~  1000  m  and  rc*  ~  2  m/s  as  representative  implies 
an  approximate  time  scale  of  about  8  min.  The  surface-layer  eddies  have 
a  time  scale  of  z/u^;  with  z  ~  10  m  and  n*  ~  0.3  m/s,  z/u^  is  about  30 
s.  Therefore  the  time  scale  of  the  convective  eddies  should  be  one  to  two 
orders  of  magnitude  larger  than  that  of  the  surface-layer  eddies.  Given  this 
disparity,  one  might  reasonably  expect  the  surface-layer  eddies  to  obtain  a 
rough  equilibrium  with  the  gusts. 

So  far,  I  have  discussed  only  large-scale  variations  in  wind  velocity.  Of 
course  other  quantities,  such  as  the  surface  heat  flux  or  temperature,  could 
also  vary  locally.  These  variations  could  even  occur  over  uniform  terrain,  as 
a  result  of  clouds  blocking  incoming  solar  radiation,  and  other  processes. 
In  this  report  I  chose  to  neglect  such  possible  variations  for  simplicity,  al¬ 
though  these  could  be  treated  in  a  manner  analogous  to  the  wind  gusts. 

2.2  Statistical  Distribution 

This  section  builds  a  general  framework  for  calculating  modifications  to 
surface-layer  statistics  caused  by  the  gusts  (large-scale  convective  eddies). 
The  basic  idea  is  to  specify  a  probability  density  function  for  the  wind- 
velocity  components  at  some  reference  height.  Local  application  of  Monin- 
Obukhov  similarity  then  allows  calculation  of  other  statistics. 

To  begin,  we  follow  Godfrey  and  Beljaars  (1991)  and  Grachev  et  al  (1998) 
by  defining  a  horizontal  gust  velocity  G  (the  contribution  of  the  large-scale 
convective  eddies  to  the  local  wind,  denoted  Wg  by  Godfrey  and  Beljaars 
and  U  by  Grachev  et  al)  through  the  relation 

IJl{z)  =  {IJ{z))  +  G{z)  .  (4) 

(Boldface  type  in  this  report  indicates  a  vector  quantity;  the  same  sym¬ 
bol  without  bolding  represents  the  magnitude  of  the  vector.)  Taking  the 
squared  magnitude  of  both  sides  of  equation  (4)  results  in 

Ul{z)  =  [{U{z))  +  GAz)f  +  Gl{z)  ,  (5) 

where  Gx  (z)  and  Gy  (z)  are  the  along-wind  and  crosswind  components  of 
the  gust  velocity,  respectively.  Averaging  both  sides,  while  keeping  in  mind 
that  (G  (z))  =  0,  results  in 

{Ul{z))  =  {U{z)f  +  {G^{z))  ,  (6) 
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where  {z))  =  (G^  (z))  +  {Gy  (z)).  Godfrey  and  Beljaars  and  Grachev 
et  al  set 


{G‘^  {Zr))  =  (3^wl  ,  (7) 

where  (3  is  an  empirical  factor,  called  the  gustiness  coefficient.  In  this  report,  I 
set  (3  =  0.95  and  Zr  =0.1  Zi  for  reasons  discussed  in  appendix  A.  For  brevity, 
I  usually  drop  the  explicit  functional  dependence  of  quantities  such  as  U  (z) 
and  G  (z)  on  z  when  they  are  to  be  evaluated  at  z  =  Zr.\  also  abbreviate 
the  local  wind  velocity  (U  {zr))^  as  \Jl,  and  (U  (zr))  as  Ur. 

Provided  that  the  gusts  have  random  orientation,  equation  (7)  implies  that 
the  variances  (G^)  and  (G^)  at  the  reference  height  must  individually  be 
(3‘^w1/2.  To  complete  a  probabilistic  model,  I  will  make  the  assumption  that 
Gx  and  Gy  are  independent,  Gaussian  random  variables.  The  approximate 
Gaussian  nature  of  turbulent  fluctuations  is  well  established  (Batchelor, 
1953,  pp  169-170),  whereas  the  independence  of  Gx  and  Gy  follows  from 
the  isotropic  nature  of  the  gusts.  The  joint  probability  density  function  (pdf) 
for  Gx  and  Gy  is  therefore 


1  /  Gl  +  Gl\ 

PG  {Gx,  Gy)  =  exp -  y  .  (8) 

7r(3^wi  y  f3^wi  J 

Given  the  pdf  for  the  gusts,  surface-layer  statistics  can  be  calculated  from 
the  general  formula 


if  {Ur,Qs,Z,Z0,Zi,g/Ts)) 


f  {Ur,  Qs,  z,  zq,  Zi,  g/Ts,  Gx,  Gy)  pg  {Gx,  Gy)  dGxdGy  , 


(9) 


where  f  {Ur,Qs,z,zo,Zi,g/Ts,Gx,Gy)  is  a  function  of  interest,  and 
(/  {Ur,  Qs,  z,  Zq,  Zi, . . . ))  its  ensemble  mean.  As  an  example,  consider  the 
mean  wind  speed  at  the  reference  height,  {Ul)-  Since 


Ul  {Ur,  Gx,  Gy)  —  \J {Ur  +  Gx)'^  -|-  Gy  , 


(10) 


we  have 


(Ul) 


pCXD  ('OQ  j -  / 

/  V  U  ^a:)^  +  Gy  exp  I  ■ 

J — oo  \ 


(3‘^wl  j 


dGxdG. 


y  ■ 


(11) 


Solution  of  this  type  of  integral  is  discussed  in  appendix  B.  The  result  is 


^  ^-uh 2 

I3w^  2 


i  +  ul 


' — 9 

u 


^2.  ful' 


(12) 


where  Ur  =  Ur/ (3w^,  and  the  In's  are  modified  Bessel  functions  of  the  first 
kind.  In  the  limit  Ur  0,  one  can  show  {Ul)  y/n / 2. 
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Of  course,  to  be  useful,  this  method  must  be  extended  to  quantities  in 
addition  to  {Ul)-  We  would  like  to  calculate,  for  example,  mean  profiles 
and  variances  in  the  presence  of  gusts.  This  is  where  the  assumption  of 
local  equilibrium  enters  into  the  problem.  The  basic  strategy  is  to  use  lo¬ 
cal  equilibrium  to  determine  a  local  friction  velocity,  u^,l,  as  a  function  of 
{Ur,  Qs,  z,  zq,  Zi,  g/Ts,Gx,Gy).  Then,  given  an  equation  for  the  quantity  of 
interest  in  terms  of  the  friction  velocity,  we  can  apply  equation  (9)  to  deter¬ 
mine  its  statistics. 

The  basic  Monin-Obukhov  similarity  equations  for  the  vertical  gradients  of 
mean  wind  and  temperature  are  the  starting  point  for  formulating  such  a 
solution.*  These  are  (e.g.,  Panofsky  and  Dutton,  1984;  Kader  and  Yaglom, 
1990) 


d{U) 

dz 

d{T) 

dz 


(13) 

(14) 


Here  T*  =  —Qs/u„  is  the  surface-layer  temperature  scale,  n  ~  0.40  is  the 
von  Karman  constant,  and  Pt  zz  0.9  is  the  turbulent  Prandtl  number.  The 
functions  and  (j)H  are  hypothesized  to  have  a  universal  dependence  on 
the  nondimensional  height  Q  =  z/L,  where  L  =  —vfTs/  {ngQs)  is  called  the 
Obukhov  length.  The  form  I  use  in  this  report  for  the  functions  (pM  and  (pn 
is  (Carl  et  al,  1973) 

(t>M,H  (C)  =  (1  —  O.M,hC)  )  (15) 


where  om  and  an  are  empirical  constants.  The  vertical  profiles,  found  by 
integrating  equations  (13)  and  (14),  are 


{U{z)) 

{T{z)) 


rt* 

K 


In - T'm  (  t  (  “T 


zq 


L 


L 


T{zh)  + 


T,Pt 


Zh 


In - 4'//  T  +  ^ 


L 


and 


Zh 


L 


(16) 

(17) 


The  quantity  zq  is  the  aerodynamic  roughness  length  (the  height  at  which 
the  velocity  vanishes),  and  Zh  is  a  reference  height  where  the  temperature 
is  measured.  The  functions  T'm  and  are  defined  by 


Jo  X 


(18) 


Performing  this  integration,  one  finds  (Grachev  et  al,  1998) 


{()  =  ?  In  _  Cardan  ^  .  (19) 


‘Throughout  this  report,  I  use  "temperature"  implicitly  to  mean  the  virtual  potential  tem¬ 
perature.  The  virtual  potential  temperature  includes  corrections  for  the  effect  of  water  vapor 
and  air  buoyancy  (Stull,  1988). 
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The  basic  premise  of  the  local  equilibrium  theory  is  that  these  equations 
apply  locally  to  the  gusts.  Hence  we  replace  them  with 


U*L 


{bx  cos  9  +  By  sin  6) - 


In - T'm 

zo 


i) 


M 


T{zh)- 


PtQs 

ku^l 


Zh 


In  -  -  ^  1  + 


Ll 


Ll 


and 


(20) 

(21) 


where  6  is  the  local  wind  direction,  is  the  local  friction  velocity  induced 
by  the  gusts,  Ll  =  (KgQs),  and  Bx  and  By  are  the  unit  vectors  in  the 

X-  and  y-directions.  Orienting  the  coordinate  system  so  that  the  mean  wind 
direction  is  along  the  x-axis,  one  has  immediately  (with  application  of  the 
constraint  (U)  = 


where 


=  i^{U{z)) 

=  0, 


and 


hM,H  {x,  y)  =  In - (x)  +  {v)  ■ 

y 


Taking  the  magnitude  of  equation  (20)  results  in 


U^L 


Zr  Zo 


UL  =  —hM[jP,-^ 
K  \Ll  Ll 

Equating  this  result  to  equation  (10)  yields 


Zr  Zo 


(22) 

(23) 

(24) 

(25) 

(26) 


This  equation  can  be  solved,  in  principle,  for  as  a  function  of  Ur,  Qs, 
Zo,  Zi,  g/Tg,  Gx,  and  Gy.  However  because  Hm  in  equation  (25)  is  a  com¬ 
plicated  nonlinear  function  of  and  is  in  turn  proportional  to  it 
is  not  possible  to  analytically  solve  equation  (25)  for  Fortunately  the 
equation  can  be  solved  numerically  for  u^,l  without  difficulty  The  simplest 
method  is  to  calculate  a  look-up  table  of  Ul  values  for  a  range  of  uniformly 
distributed  u^l  values.  Then,  given  a  new  value  of  Ul,  the  corresponding 
value  of  u,,L  can  be  found  by  interpolating  in  the  table.  Figure  2  shows  the 
dependence  of  [7^ /re*  on  u,,l/w*  for  several  values  of  zo/zi.  (The  motiva¬ 
tion  for  normalizing  the  velocities  by  re*  and  then  plotting  with  zq  j  Zi  as  a 
parameter  are  discussed  in  sect.  2.3.) 

When  temperature  statistics  are  calculated,  it  is  necessary  to  specify  Zh  and 
T  (zh)  in  equation  (21).  This  is  actually  very  problematic  in  the  context  of 
the  local  equilibrium  model.  The  most  reasonable  approach  may  initially 
seem  to  be  setting  Zh  and  T  (zh)  to  constant  (global)  values,  with  T  (zh) 
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Figure  2.  Local  wind 
speed  as  a  function  of 
local  friction  velocity. 
Four  curves  are  shown, 
corresponding  to 
different  values  of  ratio 
zojzi.  (For  Zi  =  1000  m, 
z^/zi  =  10“®  is 
representative  of  a 
grassy  surface,  whereas 
10“^  is  representative  of 
mature  crops  or  brush.) 


being  the  surface  temperature  and  Zh  the  thermal  roughness  length  of  the 
surface.  The  problem  with  this  approach  is  that  the  local  fluctuations  in 
(and  therefore  then  create  extremely  large  temperature  fluctuations 
at  z  =  Zr-  These  fluctuations  are  associated  with  turbulence  length  scales 
above  the  cutoff  for  the  local  filter;  that  is,  they  come  from  the  convective, 
Zj-scale  eddies  whose  velocity  fluctuations  are  specified  by  equation  (8).  In 
actuality,  observations  suggest  that  the  temperature  variance  cj|,  does  not 
depend  on  the  large  eddies,  even  well  above  the  surface  layer.  Therefore 
this  approach  is  not  physically  realistic.  The  root  of  the  difficulty  is  that 
Zh  and  T  (zh)  vary  in  response  to  the  gusts.  As  discussed  by  Panofsky  and 
Dutton  (1984,  p  147),  Zh  should  be  order  kh/u^,  where  kh  is  the  molecular 
coefficient  of  heat  conduction.  In  the  context  of  local  equilibrium,  then,  Zh 
is  order  kh/u^L-  When  the  wind  slackens,  the  layer  adjacent  to  the  ground 
thickens  where  molecular  heat  conduction  dominates.  A  strong  gust  ex¬ 
tends  turbulent  mixing  closer  to  the  surface,  causing  the  molecular  layer  to 
shrink. 

To  circumvent  the  problem  of  specifying  the  dynamic  response  of  Zh  and 
T  (z/i)  to  gusts,  I  will  assume  that  the  variance  of  the  locally  averaged  tem¬ 
perature  fluctuations,  {Tq)  (where  Tq  =  (T)^  —  (T)),  is  zero  at  z  =  Zr,  as  is 
consistent  with  observations.  Then  we  can  simply  apply  equation  (21)  with 
Zh  =  Zr  and  with  T  (zr)  fixed  to  the  global  value.  Results  derived  from  this 
approach  must  be  viewed  as  speculative;  fortunately  only  a  few  of  the  cal¬ 
culations  (specifically  the  temperature  gradient  and  variance)  in  this  report 
depend  on  it. 

2.3  Nondimensionalization 

Suppose  we  apply  the  methodology  in  the  previous  section  to  determine 
statistics  in  the  atmospheric  surface  layer.  Such  statistics  would  generally 
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depend  on  the  parameters  Ur,  Qs,  g/Tg,  z,  zq,  and  Zi.  (Other  parameters, 
such  as  L  and  w^,  can  be  written  as  combinations  of  this  set.)  It  would 
be  rather  clumsy,  though,  to  systematically  study  the  dependence  of  the 
statistics  on  all  five  of  these  parameters.  Therefore  we  might  ask  whether 
the  parameter  dependence  can  be  simplified  in  some  fashion.  In  this  section 
I  show  that  it  is  possible  to  reduce  the  dependence  to  just  three  basic  ratios 
involving  the  five  original  parameters,  specifically  Ur/w„,  z/zi,  and  ZQjzi. 
The  zj Zi  dependence  drops  out  if  z  is  not  part  of  the  original  parameter  set. 

As  will  become  apparent  in  section  3,  most  wind  and  temperature  statis¬ 
tics  of  interest  are  proportional  to  ensemble  averages  involving  powers  of 
the  local  friction  velocity,  along  with  cosines  and  sines  of  the  local  wind 
direction.  Therefore  we  are  faced  with  the  general  problem  of  converting  a 
function  whose  explicit  dependence  is  on  cos  d,  and  sin  d,  to  one  with 
dependence  on  Ur,  Qs,  etc,  so  that  equation  (9)  can  be  applied.  To  start, 
consider  equation  (25)  evaluated  at  z  =  Zr  and  recast  in  the  following  form: 


Ul  _  I  U^L,  (  KZrwl  HZqwI 

I 

W„  K  W„  \  ZiU-l^  ZiU'l^ 


(27) 


This  version  uses  the  relationship  =  —{ziln)  {u„l/wQ^.  Since  the  value 
0.1  Zi  is  used  for  Zr,  the  right-hand  side  is  a  function  entirely  of  u„l/w*  and 
zo/zi-  If  I  assume  that  Hm  is  an  invertible  function,  then  some  function  b 
must  exist  such  that 


U^L  _  ^  /  I/l  Zq 
W„  Zi 


Furthermore,  from  basic  trigonometric  arguments, 

Ur  -\-  Gx 


COS  6  = 


sin0  = 


^J{Ur  +  Gxf  +  Gl 

Gy 

^J{Ur  +  Gxf  +  Gl 


,  and 


(28) 


(29) 

(30) 


With  hM  given  by  equations  (19)  and  (24),  analytical  determination  of  the 
function  b  is  impossible.  The  important  point  here,  though,  is  the  existence 
of  the  function  b.  Now,  suppose  we  wish  to  know  the  average  value  of  some 
function  /  {u^l /w„  ,  9) .  Using  equations  (28),  (29),  and  (30)  to  determine  the 
dependence  of  /  on  Ur/w„,  Gxlw„,  Gy/w„,  and  z^/zi  results  in  the  proba¬ 
bility  integral 


/I— .9|  >  = 


/oo  r 

-oo  J  —  oo 


f 


arctan 


Gy  Y 
Ur  +  Gx) 


Pg  {Gx,  Gy)  dGxdGy 


(31) 
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Transforming  the  variables  of  integration  to  the  dimensionless  ratios  G'^  = 
Gxjw^  and  G'y  =  Gy/w^,  and  assuming  that  the  joint  pdf  depends  only  on 
these  ratios  (as  does  the  Gaussian  pdf  in  the  previous  section,  eq  (8)),  I  have 


arctan 


Gy  Y 

Ur  +  GxJ 


PG  {Gi,G'y)  dG'jG'y. 


(32) 


The  function  /  now  depends  entirely  on  Ur/w^,  zq/ zi,  G'^,  and  Gy.  The  de¬ 
pendence  of  (/)  on  G'^  and  G'y  is  of  course  eliminated  by  the  integration, 
leaving 


/ 


=  g 


\w^  Zi  J 


(33) 


where  g  is  the  result  of  the  integration.  This  theorem  demonstrates  that 
statistics  involving  u^l/w^  and/or  9  depend  only  on  Ur/w^  and  zojzi,  al¬ 
though  in  practice  it  may  be  impossible  to  actually  obtain  an  analytical  ex¬ 
pression  for  the  function  g{Urlw^f,zol  Zi). 

An  interesting  special  case  of  the  preceding  theorem  that  applies  to  the 
average  friction  velocity  in  the  limit  of  vanishing  mean  wind  was  proposed 
by  Businger  (1973b).  He  suggested  that 


{u^l) 


(34) 


where  gs  is  an  unknown  function.  The  validity  of  this  scaling  relationship 
was  later  supported  by  Schumann's  (1988)  analysis  based  on  a  simplified 
bulk  model  for  the  atmospheric  surface  layer.  Equation  (33)  implies  the 
same  scaling  relationship,  since  the  dependence  on  Ur/w^  disappears  in 
the  limit  of  vanishing  Ur-  In  this  report,  the  procedure  for  deriving  this 
relationship  was  much  different  from  Schumann's,  however,  being  based 
instead  on  the  statistical  nature  of  the  gusts.  Furthermore,  it  is  now  clear 
that  any  statistic  of  u^l  and/or  9  depends  solely  on  zo/zi  when  the  mean 
wind  is  zero. 

Several  valuable  corollaries  of  equation  (33)  can  be  easily  derived.  First,  let 
us  consider  statistics  of  quantities  that  additionally  depend  on  z/Ll.  Since 
z/Ll  oc  (z/zi)  arid  the  ratio  has  already  been  accounted 

for,  a  dependence  on  z/zi  is  introduced: 


/ 


U^:L 


\W^  '  Zi  Zi) 


(35) 


A  local  dependence  on  zq/Ll  would  not  add  new  dependencies  to  the  re¬ 
sult,  since  zq/L^  oc  {zoizi)  {w^/u^l)^.  Next  consider  statistics  that  depend 
onT^L  =  —Qslu^fL,  as  is  usually  the  case  for  functions  involving  tempera¬ 
ture.  Defining  a  mixed-layer  temperature  scale  as  0*  =  —Qs/w^,  we  can 
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write  Again,  no  new  dependencies  are  added  to  the  re¬ 

sult.  Hence  we  have,  most  generally. 


/  ,  f  ^  ^  _LAL/9^\-  ^ 

\  '  Ll  Ll  J  /  ^  \  Zi'  Zi 


(36) 


An  important  application  of  equation  (35)  made  later  in  this  report  is  to 
statistics  of  the  height-dependent  gust  profiles  Gx{z)  and  Gy{z).  From 
equations  (20)  and  (4),  we  find 


Gx  (z) 
Gyjz) 


1  U^L 
K 


COS  9  Hm 


■  au 

- sin  9  liM 

K 


kzqwI  \ 
>^ZqwI  \ 

Ziulh  J 


{U{z)) 


and 


1C* 


(37) 


(38) 


The  terms  in  equations  (37)  and  (38)  involving  the  function  Hm  depend  on 
u^l/w*,  z/Ll,  and  zq/zi.  Therefore,  equation  (35)  shows  that  statistics  of 
these  terms  depend  on  Ur/w^,  zjzi,  and  zq/a-  Taking  the  global  average  of 
equation  (20)  yields 


{U{z)) 

1C* 


^  /'^*L  au 
—  (  - cos  9  riM 

K  \  1C* 


/  Z  KZqW^  \  \ 

Ul’  z,nlJ/ 


(39) 


We  see  that  this  term  also  depends  on  Ur/w^,  zjzi,  and  zq/a-  Therefore 


/ 


Gx{z)  Gy{z) 


1C* 


1C* 


Vic*  '  Zi  Zi) 


(40) 


Generalized  Similarity  Theory 

Although  it  is  not  immediately  obvious,  the  nondimensionalization  theo¬ 
rems  in  the  previous  section  imply  modifications  to  the  Monin-Obukhov 
and  other  similarity  theories  for  the  surface  layer.  Before  discussing  these 
modifications,  I  will  review  several  of  the  common  similarity  theories: 

Monin-Obukhov  similarity  (Monin  and  Obukhov,  1954):  The  basic  physical  pa¬ 
rameters  hypothesized  to  be  significant  are  it*,  Qs,  b  =  g/Tg,  and  z.  This  pa¬ 
rameter  set  possesses  only  three  fundamental  physical  dimensions:  length, 
time,  and  temperature.  According  to  Buckingham's  tt  theorem  (see,  for  ex¬ 
ample,  Stull,  1988),  we  should  therefore  select  just  three  of  these  parameters 
as  key.  The  three  that  are  conventionally  selected  are  u*,  Qs,  and  z.  Quan¬ 
tities  scaled  (made  dimensionless)  by  these  parameters  must  be  functions 
of  a  single  dimensionless  ratio.  This  dimensionless  ratio  is  taken  to  be  z/ L, 
where  L  =  —u^Tg/  (ngQg)  is  the  Obukhov  length.  (The  minus  sign  and 
constant  k  in  L  are  also  matters  of  convention.) 

Local  free-convection  similarity  (Wyngaard  et  al,  1971):  The  basic  physical  pa¬ 
rameters  are  the  same  as  Monin-Obukhov  similarity,  except  that  ii*  is 
dropped  from  the  set.  The  remaining  three  parameters  incorporate  the 


three  physical  dimensions  of  length,  time,  and  temperature.  As  no  di¬ 
mensionless  combinations  can  be  formed  from  these  parameters,  quanti¬ 
ties  scaled  by  them  must  be  constant.  For  example,  quantities  with  dimen¬ 
sions  of  velocity,  when  divided  by  the  derived  scale  Uf  =  {Qszg/Tg)^^^,  are 
constant. 

Mixed-layer  similarity  (Deardorjf,  1970):  This  scheme  is  the  same  as  free- 
convection  similarity,  except  that  Zi  replaces  z.  The  velocity  scale  = 
(QsZig/Ts)^^^  therefore  replaces  Uf. 

All  three  of  the  preceding  similarity  theories  have  been  applied  to 
surface-layer  statistics.  Monin-Obukhov  similarity  is  widely  applied  to  the 
mean  profiles  for  both  mechanical  and  convective  turbulence.  Local  free- 
convection  similarity  is  applied  to  the  temperature  and  vertical  velocity 
variance  in  convective  conditions.  Mixed-layer  similarity  is  applied  to  the 
horizontal  velocity  variance  in  convective  conditions. 

Suppose  now  that  the  parameters  from  all  three  of  these  schemes  are  com¬ 
bined,  resulting  in  the  five-parameter  set  tt*,  Qs,  b,  z,  and  Zi.  Three  of  these 
parameters  (or  combinations  thereof)  can  be  taken  as  key.  Quantities  scaled 
by  the  three  key  parameters  must  depend  on  two  dimensionless  ratios.  Nu¬ 
merous  combinations  of  key  parameters  and  dimensionless  ratios  are  pos¬ 
sible.  To  generalize  Monin-Obukhov  similarity,  we  would  logically  take 
u*,  Qs,  and  z  as  the  key  parameters,  and  the  two  dimensionless  ratios  as 
z/L  and  Zi/L.  Alternatively,  we  could  formulate  a  generalized  mixed-layer 
scheme  with  Qs,  b,  and  Zi  as  the  key  parameters.  For  the  dimensionless  ra¬ 
tios,  the  only  restriction  is  that  n*  and  z  must  be  uniquely  recoverable.  So 
we  may  take  these  ratios  to  be  u„/w„  and  z/zi,  for  example. 

The  roughness  length  zq  is  not  included  in  any  of  the  three  similarity  the¬ 
ories  described  earlier.  Omission  of  zq  can  be  justified  on  the  basis  that  the 
structure  of  a  high  Reynolds  number  flow  should  be  independent  of  rough¬ 
ness  at  heights  z  >  zo  (Termekes  and  Lumley,  1972,  p  146).  However,  this 
reasoning  overlooks  possible  interactions  between  large-scale  eddies  and 
the  surface,  a  phenomenon  that  plays  a  central  role  in  the  model  described 
in  this  report.  Therefore  zq  should  generally  be  included  in  the  parameter 
set.  Forming  a  six-parameter  set  by  adding  zq  to  the  five  discussed  in  the 
preceding  paragraph,  we  conclude  that  scaled  quantities  must  depend  on 
three  dimensionless  ratios.  For  the  generalized  mixed-layer  scheme,  these 
can  be  taken  to  be  u„/w„,  z/zi,  and  zq/a- 

The  nondimensionalization  theorems  in  section  2.3  are  closely  analogous 
to  this  generalized  roughness-inclusive  mixed-layer  scheme.  The  main  dif¬ 
ference  is  that  mean  wind  speed  Ur  appears  in  place  of  u„ .  Such  a  scaling 
could  have  been  argued  by  choosing  Ur,  Qs,  b,  z,  zq,  and  Zi  as  the  initial 
six-parameter  set,  as  long  as  the  reference  height  does  not  introduce  a  new 
scale  (i.e.,  the  reference  height  must  be  proportional  to  either  zq  or  Zi).  The 
reason  for  preferring  Ur  over  rt*  in  the  present  model  is  that  it  simplifies 
the  formulation  of  the  pdf  for  the  wind  velocity  at  the  reference  height.  We 
now  see  that  the  nondimensionalization  theorems  are  consistent  with,  and 
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could  have  been  anticipated  on  the  basis  of,  a  generalized  mixed-layer  sim¬ 
ilarity  theory  that  incorporates  Ur,  z,  and  zq,  in  addition  to  the  standard 
parameters  Qs,  Zi,  and  b.  The  model  in  section  2.2  provides  a  framework 
for  calculating  the  similarity  functions. 
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3.  Effect  of  Gusts  on  Atmospheric  Statistics 


In  this  section,  I  apply  the  framework  for  calculating  gust  effects  on  at¬ 
mospheric  statistics.  The  calculations  involve  the  numerical  integration  of 
equation  (9)  with  pc  {Gx,  Gy)  given  by  equation  (8),  and  u^l  determined 
from  equation  (25).  The  numerical  integrations  were  done  on  a  square  grid 
defined  by  —12w^  <  {Gx,Gy)  <  12rc*.  The  grid  had  201  evenly  spaced 
points  in  both  the  x-  and  p-directions.  The  calculations  require  only  a  few 
seconds  each  on  a  desktop  personal  computer. 

3.1  Friction  Velocity 

Let  us  begin  by  applying  the  methodology  developed  in  section  2.2  to  the 
friction  velocity.  The  familiar  global  version  of  the  friction  velocity  is 
defined  through  the  equation 

ul  =  -  {uw)  .  (41) 

In  this  equation,  u  is  the  fluctuation  in  the  horizontal  wind  component 
aligned  with  the  mean  wind  and  w  is  the  fluctuation  in  the  vertical  com¬ 
ponent.  Similarly,  the  local  friction  velocity  is  defined  as 

uIl  =  -  ,  (42) 

where  the  prime  indicates  the  coordinate  system  aligned  with  the  local 
wind  direction  9.  By  elementary  trigonometry, 

u  =  u'  cos  9  —  v'  sin  9  ,  (43) 

where  v  is  the  horizontal  crosswind  component  of  the  wind.  Multiplying 
through  hy  —w  and  taking  the  local  average,  one  finds 

—  {uw)  =  —  (u'w') ^  cos  9  +  sin  9  .  (44) 

When  I  substitute  equation  (42)  and  set  {v'w)^^  =  0  (a  result  of  local  sym¬ 
metry  in  the  turbulence  statistics),  this  becomes 

-  {uw)^  =  uIl  cos  9  .  (45) 

Now,  by  taking  the  global  average  and  recalling  that  the  global  average 
of  a  local  average  equals  the  customary  global  average  (see  sect.  2.1),  one 
obtains  the  result 


uI  =  {uIlCos9)  .  (46) 

We  see  that  the  global  friction  velocity  is  not  the  simple  average  of  the  local 
one.  The  changing  wind  direction  plays  a  role  through  the  factor  cos  9.\{  9 
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is  randomly  distributed  in  all  directions,  the  cos  6  factor  has  zero  mean,  and 
ii*  vanishes. 

The  global  friction  velocity  is  compared  to  the  averaged  local  value,  (u^l), 
in  figure  3.  For  wind  speeds  Ur  larger  than  w^,  the  global  and  averaged  local 
values  are  nearly  the  same.  As  the  wind  speed  approaches  zero,  however, 
the  global  value  vanishes,  whereas  the  local  average  becomes  constant. 
This  asymptote  is  the  minimum  friction  velocity  recognized  by  Businger 
(1973a,b),  Schumann  (1988),  and  others.  Also  shown  on  the  figure  are  ap¬ 
proximate  results  for  (u^l)  based  on  the  average  local  wind  speed  at  the 
reference  height,  (Ul),  as  given  by  equation  (12).  I  found  these  results  by 
solving  equation  (25)  for  u^l,  except  that  I  replaced  the  local  wind  speed  by 
(Ul)-  The  resulting  value  of  will  be  designated  as  in  the  remainder 
of  this  report;  that  is. 


(Ul) 


K 


hjM 


V  -^eff  ’  Lef[ 


(47) 


where  Lefr  =  —u^^esTg/  {KgQs)-  The  figure  demonstrates  that  is  a  very 
good  approximation  to  the  actual  value  of  (u^l)  for  all  reasonable  combina¬ 
tions  of  Ur/w^  and  zq/a-  This  observation  lends  support  to  the  analyses  of 
Godfrey  and  Beljaars  (1991)  and  Grachev  et  al  (1998)  who  use  an  effective 
value  for  rt*  based  on  the  mean  wind  speed. 

Let  us  analyze  now  in  more  detail  the  vanishing  mean  wind  speed,  us¬ 
ing  approximate  forms  for  cpM  (C)  that  will  allow  us  to  obtain  an  analytical 
solution.  Simple  approximations  for  (j)M  (C)  are  available  for  the  two  limit¬ 
ing  cases  of  neutral  and  free-convection  stability.  It  may  seem  obvious,  at 
first  glance,  that  we  should  use  the  free-convection  stability  limit,  namely 
(/)M  (C)  —  (— (Garl  et  al,  1973;  Kader  and  Yaglom,  1990).  However, 
we  should  keep  in  mind  that  even  when  the  mean  wind  vanishes,  there 


Figure  3.  Global  (solid 
lines)  and  average  local 
(dashed  lines)  friction 
velocity  as  a  function  of 
vector-average  wind 
speed  at  reference  ^ 

height,  Ur  -  Dotted  lines  § 
are  an  approximation  S 
for  average  local  friction  g 
velocity  based  on  mean 
local  wind  speed,  (C/l).  "S 

lo 

E 

o 


Normalized  wind  speed,  Uflw, 
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is  still  local  shear  from  the  gusts.  Sykes  et  al  (1993)  actually  based  their 
analysis  on  locally  neutral  profiles.  So,  to  be  careful,  we  should  calculate 
the  neutral  and  free-convection  stability  limits  for  the  local  wind  profiles 
separately  and  then  compare  the  results  to  numerical  calculations  made 
with  the  full  Monin-Obukhov  similarity  equation  for  4>m  (0- 

For  the  neutral  limit,  (pM  (C)  =  1/  arid  one  has  (app  B,  eq  (B-26)) 


(w*l) 


re*  2\-D.{zr/zQ) 

Replacing  the  constants  with  their  numerical  values  results  in 


{u*l) 


=  0.34  In 


-h. 


Zi 


yiOzo 

For  the  free  convection  limit  (app  B,  eq  (B-34)), 

(u^l)  r  (5/4) 


(48) 


(49) 


^/3 


Zq) 


1/31  “1/2 


(50) 


Replacing  the  constants  with  their  numerical  values. 


=  0.42 


-2.2 

zqJ 


1-1/2 


(51) 


Generally,  Zi  is  much  larger  than  zq,  and  we  may  neglect  the  second  term 
in  the  square  brackets,  leaving 


/1s4=0,42 

w*  V  2^0  / 


(52) 


This  result  is  very  similar  to  an  approximation  derived  by  Schumann  (1988) 
(his  eq  (37)),  the  only  difference  being  the  numerical  coefficient,  which  Schu¬ 
mann  determined  as  0.52.  It  is  reassuring  that  the  approach  in  this  report, 
based  on  a  statistical  model  for  the  gusts,  yields  nearly  the  same  results  as 
Schumann's  model,  which  is  based  on  an  entirely  different  approach  in¬ 
volving  bulk  approximations. 

Figure  4  compares  a  numerical  evaluation  of  the  average  friction  veloc¬ 
ity,  based  on  local  application  of  the  full  Monin-Obukhov  equation  (15) 
for  (/m  (C)/  to  the  two  approximate  analytical  solutions,  equations  (49)  and 
(51).  Also  shown  is  rt*,efr  from  equation  (47),  which  is  observed  to  be  an 
excellent  proxy  for  (u*l)  throughout  the  range  of  roughnesses.  The  ap¬ 
proximation  that  the  profiles  under  the  gusts  locally  obey  free-convection 
similarity  produces  better  agreement  with  the  full  local  Monin-Obukhov 
evaluation  than  does  neutral  similarity.  For  unrealistically  small  roughness 
lengths  (not  shown  on  the  figure),  the  neutral  similarity  curve  does  even¬ 
tually  become  the  better  approximation.  We  conclude  that  in  conditions  of 
vanishing  wind  speed,  the  local  profiles  more  closely  resemble  their  free- 
convection  forms  than  their  neutral  ones. 
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Figure  4.  Local  average 
friction  velocity  as  a 
function  of  roughness 
length  for  vanishing 
mean  wind.  Three 
curves  are  based  on 
different  methods  for 
calculating  local  wind 
profiles  under  large 
eddies. 


3.2  Mean  Profiles 


Because  Monin-Obukhov  similarity  is  widely  used  to  model  the  mean  ver¬ 
tical  profiles  of  wind  and  temperature,  an  important  practical  matter  is 
whether  the  local-equilibrium  model  implies  any  significant,  systematic 
violations  of  global  Monin-Obukhov  similarity.  More  precisely,  when  the 
mean  gradients  are  normalized  by  u*  and  T*  as  in  equations  (13)  and  (14), 
are  they  still  universal  functions  oi  zj LI  We  might  furthermore  ask,  even  if 
the  global  application  of  Monin-Obukhov  similarity  is  tenuous,  can  empir¬ 
ically  fitted  values  for  u*  and  T*  be  devised  that  still  produce  good  agree¬ 
ment  with  atmospheric  data? 

Starting  with  the  first  question,  let  us  write  down  equations  for  the  mean 
wind  and  temperature  gradients  in  the  context  of  local  equilibrium.  Ap¬ 
plying  equations  (13)  and  (14)  locally  and  then  taking  the  global  average, 
yields 


d{U) 

dz 

d{T) 

dz 


1 

KZ 


COS  94>m 


and 


(53) 

(54) 


Next,  when  we  nondimensionalize  by  multiplying  with  nzju^  and 
—Kzu^fjPtQs,  we  have 


Kz  d  {U) 
M*  dz 
KZUi,  d  (T) 
PtQs  dz 


Plots  of  the  nondimensional  wind  and  temperature  gradients  calculated 
from  these  equations  are  shown  in  figures  5  and  6.  Four  different  values 
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Figure  5. 
Nondimensional 
wind  gradient 

{zk/u^)  {d  (U)  jdz)  as  a 

function  of  ^  =  —z/L. 


Nondimensional  height,  -z/L 


of  the  ratio  Ur/w^  are  shown:  0.3,  1,  3,  and  10.  For  a  surface  roughness 
zo/^i  =  10“*^,  these  four  values  correspond  to  —Zi/L  =  k  =  1200, 

190,  26,  and  1.6.  The  first  value  is  extremely  convective  and  probably  rare 
for  the  atmosphere.  The  second  represents  a  fair  summer  day  over  land. 

For  the  global  application  of  Monin-Obukhov  similarity  to  be  valid,  the 
curves  in  figures  5  and  6  can  depend  only  on  z/L;  that  is,  for  all  values  of 
Ur /w^f,  they  must  collapse  onto  a  single  curve.  Clearly  this  does  not  happen. 
Substantial  deviations  from  global  Monin-Obukhov  similarity  are  found 
when  Ur/w^f  <1  for  the  wind  gradient  and  Ur/w^f  <  0.3  for  the  tempera¬ 
ture  gradient.  The  dependence  of  the  nondimensional  gradients  on  Ur  / re* 
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may  explain  some  of  the  scatter  in  measurements  such  as  those  made  by 
Businger  et  al  (1971)  and  by  Hogstrom  (1988).  The  data  should  be  reexam¬ 
ined  for  dependence  on  the  inversion  height. 

Besides  Ur/w^,  the  local-equilibrium  theory  also  predicts  a  dependence  of 
the  nondimensional  gradients  on  zq/ Zi.  By  performing  numerical  evalua¬ 
tions,  I  found  this  effect  to  be  quite  insignificant,  being  most  discernable 
when  Ur/w^  <C  1. 

Next  we  turn  to  the  question  of  whether  we  can  still  fit  equations  (16)  and 
(17)  to  surface-layer  data  with  good  results,  even  when  there  is  a  break¬ 
down  in  the  global  application  of  Monin-Obukhov  similarity.  This  amounts 
to  selecting  values  for  and  T*,  possibly  different  from  their  actual  val¬ 
ues,  to  fit  the  profiles  at  some  set  of  sensor  heights.  Then  we  assess  whether 
these  fitted  values  produce  good  agreement  with  the  actual  overall  profiles. 

To  gain  some  insight  into  how  the  fitted  values  for  rt*  and  T*  can  differ  from 
their  actual  values,  consider  the  equations  for  the  ensemble  mean  profiles 
derived  from  the  local-equilibrium  theory.  The  wind  profile  follows  directly 
from  equation  (22): 

{U  (z))  =  cos  0  Hm  ■  (57) 

The  mean  temperature  profile  follows  by  averaging  equation  (21),  with 

{T{z))  =  {{T{z))l): 

(r(z))=r(zO  +  ^(r,i/rH(^,|^)^  .  (58) 

Because  hu  and  hn  are  complicated  functions  of  it  is  not  possible  to 
analytically  derive  a  general  equation  for  effective  values  of  u*  and  T*.  For 
the  idealized  situation  where  neutrally  stable  conditions  prevail  locally,  Hm 
and  hn  are  logarithmic,  and  the  two  preceding  equations  become 

{U  (z))  =  -  In  —  cos  0  )  ,  and  (59) 

K  Zo 

{T{z))  =  T{zh)-^^\n  —  {u-l)  .  (60) 

K  Zji 

Since  the  mean  profile  in  neutral  conditions,  according  to  the  global- 
equilibrium  theory,  is  {U  {z))  =  («*/«)  In  (z/zo),  the  first  of  these  equations 
implies  simply  that  (u*i,cos0)  plays  the  role  of  an  effective  value  for  tt*. 
If  the  gusts  are  relatively  strong  and  randomly  oriented,  cos  6 )  can  be 
much  smaller  than  Similarly,  —Qs  {u~l)  plays  the  role  of  an  effective 
temperature  scale  in  place  of  T*  =  —Qg/u^  in  the  global  theory. 

Turning  now  to  the  other  limiting  approximation,  where  free  convection 
prevails  locally,  one  can  show  (app  B)  that  {U  (z))  oc  cos  6  ).  Since  the 
global-equilibrium  model  predicts  {U  (z))  oc  u^,  and  ul  =  (ul^cosO)  (eq 
(46)),  the  effective  friction  velocity  in  free  convection  is  the  same  as  the 
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actual  friction  velocity.  In  regards  to  the  temperature  profile,  it  becomes 
independent  of  u*  in  free  convection  (see,  for  example  (Kader  and  Yaglom, 
1990)).  Therefore  in  this  extreme,  there  is  no  local  variability  in  the  tempera¬ 
ture  profiles,  and  T*  =  —  Qs/u*  can  be  used  as  an  effective  temperature 
scale. 

To  elucidate  these  ideas,  let  us  consider  an  example  atmospheric  boundary 
layer  with  =  3  m/s,  zq  =  10“®  m,  and  Zi  =  1000  m  {Qg  =  0.81  Ks/m). 
The  wind  speed  Ur  is  varied.  For  each  value  of  Ur,  fitted  values  of  and  T* 
are  determined  by  matching  equations  (13)  and  (14)  to  the  wind  speed  and 
temperature  at  10-m  height.  Figure  7  shows  the  resulting  ratio  of  the  fitted 
value  of  rt*  (n*,fit)  to  the  actual  value.  For  large  Urjw^,  the  fitted  friction 
velocity  is  the  same  as  the  actual  value.  As  Ur/w„  becomes  less  than  about 
10/  becomes  less  than  the  actual  value,  as  would  be  anticipated  from 
the  preceding  discussion.  Surprisingly,  as  Ur/w„  0,  the  ratio  approaches 
a  constant  value  smaller  than  one.  Based  on  the  preceeding  discussion,  the 
asymptote  should  equal  one.  The  reason  behind  this  apparent  discrepency 
is  that  the  free-convection  limit  is  never  truly  attained.  The  gusts  prevent 
the  turbulence  from  being  in  such  a  state  locally. 

Figures  8  and  9  illustrate  agreement  between  wind  and  temperature 
profiles  from  the  local-equilibrium  model  to  those  based  on  fitted  val¬ 
ues  for  u„  and  T*.  Plotted  are  the  ratios  [f/gt  {z)  —  {U  (z))]  /  {U  (z))  and 
[Tfit  {z)  —  {T  (z))]  /  (T  (z))  for  z  =  2  m  and  z  =  100  m,  where  f/gt 
and  Tgt  (z)  are  the  profiles  calculated  by  fitting  tt*  and  T*  to  measure¬ 
ments  at  z  =10  m.  Also  shown  are  [t/giob  (z)  —  {U  (z))]  /  {U  (z))  and 
[Tgiob  (z)  —  {T  (z))]  /  (T  (z))  for  z  =  2  m  and  z  =  100  m,  where  t/giob  (z)  and 
^giob  (z)  are  the  profiles  calculated  from  the  global-equilibrium  model.  The 
profiles  calculated  from  the  fitted  values  are  always  in  good  agreement 
(relative  error  less  than  about  4%)  with  the  local-equilibrium  model,  even 


Figure  7.  Ratio  of 
friction  velocity 
determined  by  fitting 
measurements  at  10-m 
height  to  actual  value 
determined  from 
local -equilibrium 
model. 
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Figure  8.  Relative  error 
(in  %)  between  wind 
speed  profile  in 
local -equilibrium  model 
and  a  profile 
defermined  by  fitting  u* 
and  r*  fo  measuremenfs 
af  10-m  heighf.  Errors 
are  shown  for  2-  and 
100-m  heighfs.  Also 
shown  are  relafive 
errors  for  global 
Monin-Obukhov 
similarify  based  on 
acfual  values  of  m* 
and  T*. 
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Figure  9.  Same  as 
figure  8,  excepf  for 
femperafure  profile. 


Mean  wind  speed,  U/zvt 


for  very  small  Ur/w^.  The  global-equilibrium  profiles  for  the  wind  speed 
are  in  poor  agreement  with  the  local-equilibrium  model,  although  agree¬ 
ment  for  the  temperature  profiles  is  quite  good.  The  reason  for  this  behav¬ 
ior  was  mentioned  earlier:  the  local  temperature  profiles  are  independent 
of  n*  in  highly  convective  conditions,  and  therefore  the  global  and  local 
models  match  in  this  case.  In  interpreting  the  results  for  the  wind  profiles, 
one  should  keep  in  mind  that  values  for  Ur/w^  <1  are  unusual  for  the  at¬ 
mosphere.  Even  for  Ur/w^  =  1,  though,  the  global-equilibrium  model  can 
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be  off  by  more  than  10  percent.  Application  of  global  similarity  theory  to 
the  wind  profiles  in  highly  convective  atmospheric  conditions  is  therefore 
very  tenuous,  although  when  fitted  values  for  n*  and  T*  are  used  in  place 
of  the  actual  ones,  the  theory  can  be  made  to  work  fairly  well. 


3.3  Temperature  Variance 


In  this  section,  I  explore  the  implications  of  the  local-equilibrium  model  for 
the  temperature  variance.  The  temperature  variance  is  defined  as 

a^r{z)  =  {[T{z)-{T{z))f)  .  (61) 

Previous  studies,  based  on  a  global  similarity  approach,  have  found  the 
surface-layer  temperature  variance  in  neutral  conditions  to  be  about  4.0T^^ 
(e.g.,  Stull's  (1988)  equation  9.5.3k).  In  freely  convective  conditions,  the  pre¬ 
vious  studies  suggest  =  0.90T^^  {—zjL)~‘^^^  (Kaimal  et  al,  1976;  Caughey 
and  Palmer,  1979).  The  following  equation  smoothly  interpolates  between 
these  two  limits: 


(T4 

^=40 

J<2  ^ 


-1 


(62) 


Note  that  the  temperature  fluctuations  are  entirely  "small  scale";  that  is, 
since  neither  tc*  nor  Zi  appear  in  the  previous  equation,  the  large-scale  con¬ 
vective  eddies  have  not  previously  been  observed  to  affect  the  temperature 
variance. 


In  the  context  of  local  equilibrium,  the  large  convective  gusts  can  modify 
the  temperature  variance  in  two  ways.  First,  they  can  directly  vary  the  tem¬ 
perature  through  local  adjustments  to  the  vertical  profile.  Second,  they  can 
indirectly  affect  the  temperature  variance  by  locally  modulating  through 
changes  in  u^l.  Let  us  decompose  the  temperature  as  follows: 

T  =  {T)+TG  +  t,  (63) 


where  Tq  is  the  change  induced  directly  by  the  gusts  and  t  represents  the 
small-scale  fluctuations  obeying  equation  (62),  except  with  local  scales  re¬ 
placing  the  global  ones.  The  variance  is  therefore 

4  =  ((Tg  +  tf)  =  (4)  +  2  {Tot)  +  (f)  .  (64) 

Let  us  decompose  t  as 

t  =  CFT,Li  ,  (65) 

where  ^  is  a  random  variable  having  zero  mean  and  unit  variance.  The  fac¬ 
tor  iTr,L  scales  the  variance  locally  in  response  to  the  gusts.  Having  factored 
out  crr,L,  I  can  now  consider  the  random  variable  ^  a  process  independent 
from  the  gusts.  As  a  result,  the  direct  and  indirect  contributions  separate 
from  each  other: 


4  —  (Tg)  -I-  2  {Tg<jt,l)  iO  +  —  (Tg)  +  ■  (66) 
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Let  us  first  consider  the  direct  contribution,  By  definition,  Tq  = 

(r  {z))]^  —  {T  (z)).  From  equation  (21)  (with  Zr  appearing  in  place  of  Zh, 
as  discussed  at  the  end  of  sect.  2.2), 


rp  Pt 

=  — 

K 


(67) 


and  therefore 


Substituting  T^l  =  —QsIu^l,  and  normalizing  with  0*  =  Qs/w^,  we  have 


Since  Zr/L^  is  proportional  to  (vj^/u^l)^,  the  right-hand  side  is  a  function 
of  and  z/L,  and  we  know  from  the  theorem  (35)  that  {Tq)  /OI  must 

depend  on  the  dimensionless  ratios  Urlw^,  z/zi,  and  zq/  Zi. 

With  respect  to  the  indirect  contribution,  the  local  temperature  variance  is 
the  same  as  equation  (62),  except  that  replaces  T*  and  replaces  L. 
Making  these  replacements  and  normalizing  the  local  friction  velocity  by 
re*,  we  have 


^T,L 


=  4.0 


1+4A 


(70) 


Similarly  to  the  direct  contribution,  we  know  from  equation  (35)  that 
{^T  l)I^*  naust  depend  on  the  dimensionless  ratios  Ur/w^,  zjzi,  and  z^j Zi. 

The  local-equilibrium  model  prediction  for  Ur/w^  =  0.3  is  broken  down 
and  compared  to  global  similarity  in  figure  10.  The  direct  effect  of  the  gusts 
(local  modification  of  the  vertical  temperature  profile)  is  insignificant.  The 
variance  is  almost  entirely  locally  generated.  But  more  importantly,  the  lo¬ 
cal  and  global  similarity  models  are  nearly  identical,  even  though  this  case 
is  highly  convective  by  atmospheric  standards.  This  behavior  can  be  ex¬ 
plained  by  the  fact  that  the  convective  contribution  to  the  variance,  0.90 

T‘^l{~zI =  0.90  {TqQII ngz)^^'^ ,  is  independent  of  u^l-  Therefore 
the  gusts  have  a  vanishing  local  effect  in  highly  convective  conditions.  I  ex¬ 
amined  several  other  combinations  of  Ur  / re*  and  zq  / Zi  and  concluded  that 
the  gusts  cause  no  significant  modifications  to  the  temperature  variance 
for  ordinary  atmospheric  conditions.  The  local-equilibrium  model  is  there¬ 
fore  consistent  with  atmospheric  observations,  which  have  so  far  shown 
no  discernable  dependence  of  the  temperature  variance  on  the  large-scale 
convective  eddies. 


Figure  10.  Temperature 
variance  predictions  for 
Ur/wi,  =0.3.  Solid 
curve  is  total  prediction 
from  local-equilibrium 
model.  Dashed  curve  is 
confribufion  fo 
local -equilibrium  model 
caused  direcfly  by 
gusfs.  Doffed  curve  is 
fofal  prediction  from 
global-equilibrium 
model. 


Temperature  variance,  <o|(z))/e? 


3.4  Vertical  Velocity  Variance 

Previous  experimental  studies  have  shown  that  the  vertical  velocity  vari¬ 
ance,  like  the  temperature  variance,  is  concentrated  at  small  scales.  Panof- 
sky  et  al  (1977)  found  the  vertical  velocity  variance  in  the  surface  layer  to 
be  well  represented  by 


=  1.6  +  2.9  -- 
ui  V  L/ 


Z\2/3 


or,  equivalently,  by 


^  =  1.6  — 


wt 


u* 


w* 


2/3 


(71) 


(72) 


Unlike  the  mean  temperature  profile,  the  mean  vertical  velocity  profile  ex¬ 
hibits  no  height  dependence:  it  is  essentially  zero  at  all  heights  in  the  sur¬ 
face  layer.  Therefore,  in  developing  a  local-equilibrium  model,  I  should  set 
the  large-scale  (direct)  contribution  to  the  vertical  velocity  to  zero.  We  have 
simply  cj^  =  where  ^  is  given  by  equation  (72),  with  replac¬ 

ing  tt*: 


4  =  1.6%^ +  1.6 

wi  w% 


(73) 


Figures  11  and  12  show  the  vertical  velocity  variance  calculated  from  the 
global-  and  local-equilibrium  models.  The  figures  are  for  roughness  lengths 
of  zq  =  10~^Zi  and  zq  =  10~^Zi,  respectively.  The  ratio  of  the  local  pre¬ 
diction  (eq  (73))  to  the  global  prediction  (eq  (72))  is  plotted.  The  local  and 
global  predictions  start  to  differ  substantially  when  tU/rc*  <  1,  particularly 
for  the  rough  surface.  In  fact,  the  variance  nearly  doubles  near  the  ground 
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Figure  11.  Ratio  of 
vertical  velocity 
variance  (crj,)  predicted 
by  local-equilibrium 
model  to  global 
similarity  prediction, 
for  a  smooth  surface 
(zo/2i  =  10-5). 
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Figure  12.  Same  as 
figure  11,  but  for  a 
rough  surface 
{zo/Zi  =  10“5). 
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for  the  rough  surface  when  Ur/w^  =  0.3.  The  cause  is  mechanical  turbu¬ 
lence  generated  locally  by  gusts. 

Let  us  analyze  the  behavior  of  the  vertical  velocity  variance  for  small  Ur  / tc* 
in  more  detail.  Recall  from  section  3.1  that,  in  analyzing  statistics  of  u^l, 
the  free-convection  approximation  for  the  profiles  was  found  to  be  more 
satisfactory  than  the  neutral  approximation.  Applying  the  free-convection 
approximation  to  the  local  wind  profiles,  we  have  from  equation  (B-34) 

wl  6 

When  I  substitute  numerical  values  for  the  constants  and  make  the  addi¬ 
tional  approximation  zq  <C  Zr,  the  previous  equation  becomes 


(75) 


Hence,  from  equation  (73), 


^  =  0.31^^^  -h  1.6 
wi  wi 


(76) 


The  following  approximate  expression  interpolates  smoothly  between  the 
mechanical  and  free-convection  limits: 


al,  =  l.Qul  + 


0.31  ( 

ZiJ 


1/3  /  N  2/3' 

+  1.6  (- 
Zi 


(77) 


The  only  difference  between  this  equation  and  the  one  from  Panofsky  et  al 
(1977)  (eq  72)  is  the  extra  term,  0.31  w‘1  {zQ/zi)^^^.  This  formula  is  plotted 
and  compared  to  the  local  equilibrium  prediction  in  figure  13.  We  see  that 
the  approximate  formula  works  best  for  /tc*  <C  1  but  is  still  satisfactory 
for  moderate  Ur/w^. 


3.5  Horizontal  Velocity  Variance 

Analogously  to  the  temperature,  the  velocity  component  in  the  along-wind 
(x-)  direction  can  be  decomposed  as  follows: 

U  =  {U)  +  G^  +  u,  (78) 

where  u  contains  the  small-scale  fluctuations  filtered  out  by  the  local  aver¬ 
aging  operation.  Therefore  we  have  for  the  variance 

=  ((H  -  {U)f)  =  (Gl)  +  2  {G^u)  +  (^2)  .  (79) 

In  principle,  we  could  devise  a  formula  for  the  horizontal  velocity  variance 
similar  to  equation  (62)  for  the  temperature  variance,  interpolating  between 
the  neutral  and  free-convection  limits.  However,  the  free-convection  limit 
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Figure  13.  Comparison 
of  approximate  formula, 
equation  (77)  (dashed 
lines),  to 

local -equilibrium  model 
prediction  (solid  lines). 
Calculations  are  for  a 
rough  surface 
(zo/zi  =  10“3). 


Vertical  velocity  variance,  g^{z)Iw^ 


for  the  horizontal  velocity  variance  has  not  yet  been  conclusively  observed 
in  the  atmosphere.  Apparently,  the  contributions  to  the  variance  from  me¬ 
chanical  turbulence  and  the  large-scale  convective  eddies  nearly  always 
dominate  any  contributions  from  the  smaller  (z-)  scale  convective  eddies. 
(See,  for  example,  the  discussion  in  Kader  and  Yaglom  (1990,  pp  651-652).) 
Therefore  only  the  contribution  from  purely  mechanical  turbulence  (the 
neutral  limit)  need  be  included  in  the  formulation  for  u.  Locally,  these  fluc¬ 
tuations  scale  in  proportion  to 

Variances  for  shear  turbulence  also  depend  on  the  direction  relative  to  the 
mean  wind.  Representative  values  for  the  along-wind  and  crosswind  stan¬ 
dard  deviations  are  2.4  u*  and  1.9  u*,  respectively  (Panofsky  and  Dutton, 
1984).  In  the  context  of  local  equilibrium,  we  therefore  model  u  as 

tt  =  (2.4cos0 -F  1.9sin0)  ,  (80) 

where  ^  is  a  zero-mean,  unit-variance  random  variable  as  in  the  previous 
section.  Based  on  the  assumption  that  ^  is  independent  of  the  gusts,  the 
cross  term  in  equation  (79)  vanishes,  leaving 

(tI  {z)  =  {Gl{z)) +  h.S{ulL  cos^  e) +  ?,.Q{ulLs\n‘^e)  .  (81) 

A  similar  treatment  for  the  crosswind  variance  results  in 

al  (z)  =  {Gl  (z))  -F  5.8  sin  +  3.6  {ul^cos  ^6*)  .  (82) 

From  equations  (35)  and  (40),  we  know  that  (Gl{z))  and  {Gy{z))  in 
equations  (81)  and  (82),  after  nondimensionalization  with  w^,  depend 
on  the  ratios  Ur/w^,  zjzi,  and  zo/zj.  The  quantities  6)  /w^  and 

/w"^  are  independent  of  height  and  are  therefore  functions  of 
Ur/w^  and  zq/z*.  We  conclude  that  cr^  (z)  /wl  and  (z)  /rcj  depend  over¬ 
all  on  Ur/w^,  zj Zi,  and  zq/zx 
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If  we  were  to  neglect  the  height  dependence  of  (z))  and  (Gy  (z))  and 
set  and  9=0,  the  previous  expressions  would  reduce  to  = 

OAbwl  +  5.8u^  and  =  OAbwl  +  3.6u^.  This  simple  addition  of  variances 
from  large-scale  convective  gusts  and  mechanical  turbulence  has  been  pre¬ 
viously  suggested  by  Panofsky  et  al  (1977),  Hojstrup  (1982),  Wilson  and 
Thomson  (1994),  and  Peltier  et  al  (1996).  But  equations  (81)  and  (82)  in¬ 
corporate  additional  phenomena  not  present  in  the  previous  formulations. 
Most  significantly,  they  include  height  dependence  of  the  gust  structure, 
and  the  conversion  of  energy  from  the  gusts  to  local  mechanical  turbulence. 

To  better  understand  the  behavior  of  equations  (81)  and  (82),  let  us  de¬ 
termine  explicit  results  for  the  idealized  cases  of  local  neutral  and  free- 
convection  profiles  underneath  the  gusts.  For  the  neutral  limit,  equation 


(20)  yields  the  two  component  equations 

11  T  Z 

— ^In  —  cos  9  =  {U{z))  +  Gx{z)  and  (83) 

K  Zo 

In  —  sin  0  =  Gy  (z)  .  (84) 

K  Zo 

Squaring  both  equations  and  taking  expected  values,  we  have 

(G^(2:))  =  ^  In^  —  (ii^^cos^  0)  —  (1/(2:))^  and  (85) 

K  Zo 

{Gl{z))  =  —  {uIl  sin^  9)  ,  (86) 

K  Zq 

where 

(U  (z)}  =  -  In  —  (u^L  cos  9}  .  (87) 

K  Zo 


The  various  moments  involving  and  9  are  evaluated  in  appendix  B  as 
equations  (B-20),  (B-22),  and  (B-23).  Incorporating  those  results,  we  have 


iz)  = 


5.8k‘^U^  ^ 

\n^{Zr/zo)  2W‘{Zr/zo) 


[In^  (2/20) -6  9.4k^'  . 


A  similar  treatment  for  the  crosswind  variance  results  in  the  same  formula 
but  with  3.6  replacing  5.8  in  the  first  term.  In  equation  (88),  the  first  term  on 
the  right  arises  from  shear  production  by  the  mean  wind.  The  first  and  sec¬ 
ond  terms  in  the  square  brackets  represent  the  direct  effect  of  the  gusts  and 
local  shear  production,  respectively.  Note  that  the  local  shear  production  is 
height  independent. 

Now  let  us  reanalyze  the  problem  in  the  free-convection  limit.  If  the  mean 
wind  speed  is  negligible,  cos^  9)  and  sin^  9)  must  both  equal 
{uIl)  /2.  leaving 

al,{z)=0A{G\z))+A.7{ulL)  ■  (89) 

Away  from  the  influence  of  the  ground,  one  would  have  ^  ss 
0.5  (G^  (zr ) )  =  0.45rc^.  Very  near  the  ground,  however,  surface  drag  dimin¬ 
ishes  the  gust  velocity,  so  that  (G^  {z) )  must  become  small.  Shear  forces  can 
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counteract  this  decrease  in  the  variance  by  locally  enhancing  (uIl)-  In  the 
limit  of  large  —z/Ll  and  —zq/Lj^,  one  has  from  equations  (13)  and  (15) 


so  that 


Ul{z) 


G  {z)  ~ 


Hl 


{G\z)) 


(90) 


(91) 


Therefore,  to  evaluate  equation  (89),  we  need  and  These  can 

be  found  in  appendix  B  as  equation  (B-34).  One  subsequently  finds  that 


(z) 


t 

2 


1  -  {zp/zf^^ 
1  -  {zo/Zrf^'^ 


Substituting  numerical  values  for  the  constants,  we  have 


CJ, 


U^V 


{z) 


=  0.45 


wt 


1  -  {zp/zf^^ 

1  -  {zQlZr)^/'^ 


f  Zr)\ 

+  0.89  — 

\ZiJ 


1/3 


(92) 


(93) 


The  first  term  on  the  right  is  the  direct  effect  of  the  large-scale  gusts, 
whereas  the  second  term  represents  local  shear  production.  These  terms 
are  analogous  to  the  second  and  third  in  equation  (88).  Equation  (93)  dif¬ 
fers  from  the  usual  one  for  the  horizontal  velocity  variance  in  a  convective 
boundary  layer,  namely  cr^  y  (z)  =  0.45w^.  These  differences  are  significant 
near  the  ground,  particularly  over  rough  surfaces.  If  zq  =  10~^ Zi,  the  lo¬ 
cal  shear  production  term  equals  0.089  w^,  which  amounts  to  a  20  percent 
increase  in  the  variance  at  z  =  Zy. 

Let  us  now  consider  some  example  calculations.  Figures  14  and  15  show  the 
ratio  of  the  along-wind  velocity  variance  predicted  by  the  local  application 
of  Monin-Obukhov  similarity  to  the  global  prediction,  =  t)Ahw‘l  +  5.8tt^. 
The  former  plot  is  for  a  relatively  smooth  surface  {zq/zi  =  10“^),  whereas 
the  latter  is  for  a  rough  surface  {zo/zi  =  10“^).  The  local  theory  predicts  a 
substantial  decrease  of  the  velocity  variance  near  the  ground,  particularly 
over  a  rough  surface  when  Ur/w^  <  1.  The  primary  cause  is  the  previously 
mentioned  slowdown  in  the  velocity  of  the  gusts  when  they  interact  with 
the  surface. 

Figures  16  and  17  show  the  ratio  of  the  along-wind  to  the  crosswind  ve¬ 
locity  variance  (u^/cr^),  for  the  smooth-surface  case  {z^/zi  =  10“®).  The 
former  plot  is  the  ratio  as  a  function  of  height,  whereas  the  latter  is  the  ra¬ 
tio  as  a  function  of  Urjw^.  In  highly  convective  conditions  <C  1), 

the  turbulence  structure  is  nearly  isotropic  and  this  ratio  converges  to  1.  In 
high-wind  conditions  {Urlw^f  »  1),  this  ratio  is  dominated  by  shear  turbu¬ 
lence  and  is  roughly  equal  to  5.8/3.6  =  1.6.  Furthermore,  the  ratio  is  larger 
near  the  ground,  indicating  a  more  pronounced  anisotropic  structure  there. 
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10-1 


Figure  14.  Ratio  of 
along-wind  velocity 
variance  (cr^)  predicted 
by  local-equilibrium 
model  to  global 
similarity  formula  0.45 
+  5.8m*,  for  a  smooth 
surface  {zo/zi  =  10“®). 
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Figure  15.  Same  as 
figure  14,  but  for  a 
rough  surface 
(zo/zi  =  10“®). 
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Figure  16.  Ratio  of 
along-wind  to 
crosswind  velocity 
variance  in 

local -equilibrium 
model,  as  a  function  of 
height,  for  a  smooth 
surface  {z^lzi  =  10“®). 


Figure  17.  Ratio  of 
along-wind  to 
crosswind  velocity 
variance  in 

local -equilibrium 
model,  smooth  surface 
case  (zo/^i  =  10“®). 
Ratio  versus  wind 
speed  is  plotted. 


Normalized  wind  speed,  Uylw, 


At  greater  altitudes,  the  large  gusts  play  a  more  prominent  role,  thereby 
diminishing  the  anisotropy. 

Figure  17  suggests  that  near  the  ground,  for  moderate  values  of  Ur/w^, 
the  crosswind  variance  may  actually  exceed  the  along-wind  variance.  This 
behavior  is  caused  primarily  by  the  action  of  gusts  whose  orientation  runs 
counter  to  the  prevailing  wind  direction.  In  the  presence  of  such  gusts,  the 
local  conditions  can  very  closely  approximate  free  convection.  Gusts  in  the 
prevailing  crosswind  direction  cannot  cause  local  slackening  of  the  wind  to 
such  a  degree. 
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As  illustrated  by  figures  14  through  17,  the  local-equilibrium  model  pre¬ 
dicts  behavior  of  the  surface-layer  horizontal  velocity  variances,  which  is 
more  complex  than  previously  thought.  It  would  be  desirable  to  have  a 
simple  equation  that  approximates  the  full  behavior  of  the  local  equilib¬ 
rium  model.  We  discovered  in  deriving  equations  (88)  and  (93)  that  there 
are  three  contributions  to  the  variance:  global  shear  production  by  the  mean 
wind,  local  shear  production  from  the  gusts,  and  the  direct  contribution 
from  the  gusts.  For  the  along-wind  variance,  the  global  shear  production 
contribution  can  be  approximated  as  5.8  u^.  This  term  would  naturally  di¬ 
minish  with  decreasing  mean  wind  speed.  Local  shear  from  the  gusts  is  im¬ 
portant  primarily  when  Ur/w^  <C  1.  Since  we  found  in  section  3.1  that  the 
free-convection  approximation  works  reasonably  well  in  such  situations, 
we  can  approximate  this  term  with  its  value  from  equation  (93),  namely 
0.89  wl  {zo/zi)^^^.  For  the  direct  contribution  from  the  gusts,  we  might  try 
using  {Ul  {z))  to  approximate  the  general  shape  of  the  profile: 


{UL{z)f 

2  {ULiZr)f’ 


(94) 


One  can  show,  in  fact,  that  this  equation  is  exact  for  both  the  local  neutral 
and  local  free-convection  profiles.  The  approximation  can  be  applied  by 
calculating  {Ul  (zr))  from  equation  (12),  solving  equation  (47)  for  u^^eS,  and 
then  using  this  value  to  estimate  {Ul  (z)).  Putting  these  results  together 
yields: 


al  {z)  =  5.8u^  -h  { 


0.45 


hM 
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(95) 


An  analysis  of  the  crosswind  variance  results  in  the  same  formula  but  with 
3.6  replacing  5.8  in  the  first  term.  This  approximation  is  compared  to  the  full 
local  equilibrium  model  in  figure  18  for  Ur/w^  =  0.3  and  zq/ Zi  =  10“®.  Also 
plotted  are  curves  based  on  local  neutral  (eq  (88))  and  local  free-convection 
(eq  (88))  profiles  under  the  gusts.  The  approximate  formula,  equation  (95), 
is  quite  close  to  the  full  Monin-Obukhov  calculation.  I  have  analyzed  other 
combinations  of  Ur  / and  zq  / Zi  with  similar  satisfactory  results.  The  ef¬ 
fective  value  for  u^l  works  quite  well  as  a  shortcut  to  the  full  calculation 
with  the  pdf. 
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Figure  18.  Local 
equilibrium  predictions 
for  along -wind  velocity 
variance  using  several 
methods  of  calculation. 
For  all  curves,  C/r  /  w*  = 
0.3  and  zo/^i  =  10“®. 
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4.  Summary 


The  local-equilibrium  model  for  gusts  developed  in  this  report  suggests 
significant  violations  of  widely  used  similarity  theories,  particularly  when 
the  atmosphere  is  very  convective.  For  example,  the  inversion  height  (zi) 
appears  to  affect  scaling  of  the  mean  wind  and  temperature  gradients.  The 
primary  physical  reason  for  this  behavior  is  that  the  gusts  create  significant 
local  modifications  to  the  friction  velocity.  On  the  other  hand,  the  rough¬ 
ness  length  appears  unimportant  for  the  gradients,  as  had  been  previously 
concluded. 

Regarding  the  variances,  the  vertical  velocity  variance  in  convective  con¬ 
ditions  is  predicted  to  have  significant  dependence  on  both  the  inversion 
height  and  roughness.  Near  the  ground,  the  vertical  velocity  variance  is  en¬ 
hanced  by  local  shear  in  the  gusts.  The  horizontal  velocity  variances  have  a 
strong  dependence  on  height  and  roughness,  even  for  very  moderate  con¬ 
vection,  that  has  not  been  previously  recognized.  The  physical  reason  is 
that  the  wind  speed  in  the  gusts  must  diminish  very  near  the  surface,  ob¬ 
taining  a  roughly  logarithmic  height  dependence.  Of  the  variances  consid¬ 
ered  in  this  report,  only  the  temperature  variance  appears  to  be  nearly  un¬ 
affected  by  the  gusts. 

The  results  of  this  report  are  in  qualitative  agreement  with  the  large-eddy 
simulations  (LES)  of  Sykes  et  al  (1993)  and  Kharma  and  Brasseur  (1997).  The 
former  authors  demonstrated  a  dependence  of  boundary-layer  structure 
on  roughness  in  free  convection.  The  latter  found  that  velocity  statistics  de¬ 
viate  significantly  from  Monin-Obukhov  similarity,  whereas  temperature 
statistics  do  not.  But  these  LES  studies  cannot  be  viewed  as  providing  con¬ 
clusive  evidence  in  support  of  the  local-equilibrium  model:  the  boundary 
condition  in  LES  is  essentially  the  same  as  the  one  in  the  local-equilibrium 
model,  except  that  it  is  applied  between  the  lowermost  grid  level  and 
the  surface,  instead  of  Zr  =  0.1  and  the  surface.  Also,  LES  has  well- 
documented  shortcomings  near  the  surface,  because  of  its  finite  resolution. 
Therefore  convincing  verification  of  the  predictions  of  the  local-equilibrium 
model  will  ultimately  require  real  atmospheric  data. 
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Appendix  A.  Reference  Height  and  Gustiness  Coefficient 


An  important  modeling  issue  is  what  value  should  be  used  for  the  reference 
height.  Grachev  et  al  (1998)  used  Zr  =  10  m,  because  this  height  is  a  stan¬ 
dard  one  for  tower-based  wind  measurements.  On  the  other  hand,  because 
10  m  is  usually  well  within  the  surface  layer,  the  structure  of  the  gusts  could 
be  significantly  disrupted  by  the  ground  there.  Schumann  (1988)  based  his 
model  on  a  surface-layer  height  of  Zr  =  O.lzi.  Sykes  et  al  (1993)  took  the 
additional  step  of  deriving  an  analytical  expression  for  Zr  {zg  in  their  no¬ 
tation).  However,  Grachev  et  al  (1997)  subsequently  demonstrated  that  the 
analysis  of  Sykes  et  al,  by  assuming  that  the  local  profiles  under  the  large 
eddies  are  logarithmic,  leads  to  implausible  results. 

I  follow  Schumarm  in  this  report  by  simply  setting  Zr  =  O.lzi.  The  main  mo¬ 
tivation  for  using  this  value  instead  of  10  m  is  that  by  placing  the  reference 
height  in  the  vicinity  of  the  top  of  the  surface  layer,  the  probability  distribu¬ 
tion  more  reasonably  represents  the  forcing  from  the  gusts  rather  than  the 
response  to  them.  Then  the  isotropic  distribution,  equation  (8),  is  justifiable. 
Admittedly,  in  conditions  with  high  mean  wind  shear  (large  Ur/wf),  the 
assumption  of  isotropy  will  break  down  even  at  Zr  =  O.lzi.  However,  the 
effect  of  the  gusts  on  surface-layer  statistics  is  negligible  in  such  situations, 
since  shear  instability,  rather  than  the  large-scale  convective  eddies,  will 
dominate  the  surface  layer.  The  prescribed  pdf  for  the  gusts  is  at  issue  only 
for  small  and  moderate  wind  speeds,  for  which  the  assumption  of  isotropy 
Zr  =  O.lzj  is  quite  reasonable. 

The  other  model  constant  is  the  gustiness  coefficient  (5.  Previous  authors 
have  suggested  several  possible  values.  Based  on  their  analysis  of  marine 
data  at  a  10-m  height,  Fairall  et  al  (1996)  concluded  that  (3  ~  1.25.  Schu¬ 
mann  (1988)  determined  that  (3  =  0.7  ±  0.1  from  his  model  analysis.  Bel- 
jaars  (1994)  deduced  that  /3  ~  1.2  from  the  numerical  simulations  of  Sykes 
et  al(1993). 

A  value  for  j3  can  also  be  inferred  from  other  sources  based  on  the  reason¬ 
ing  that  the  horizontal  velocity  variance  for  z  >  O.lz^  is  determined  primar¬ 
ily  by  the  large-scale  convective  eddies.  From  measurements  made  in  the 
lower  boundary  layer  by  aircraft,  Gaughey  and  Palmer  (1979)  estimated  the 
horizontal  velocity  variances  as  O.SSrc^.  Since  the  variances  {Gf)  and  (G^) 
are  equal  to  fP'vfl/'l,  this  implies  that  (3  ~  0.84.  Similarly,  Deardorff  and 
Willis  (1985)  found  ^  (z)  =  OAbw"^  in  their  laboratory  experiment,  which 
implies  that  (3  ~  0.95. 

For  the  calculations  in  this  report,  I  set  /3  to  0.95,  since  this  value  falls 
roughly  in  the  middle  of  the  range  of  previous  research  results  reviewed 
here. 
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Appendix  B.  Moments  of  Friction  Velocity  and  Wind  Direction 


This  appendix  derives  equations  for  various  moments  of  the  friction  veloc¬ 
ity  and  wind  direction  that  appear  in  the  main  text.  The  moments  have  the 
general  form  cos™  0  sin"^  0).  Unfortunately,  as  discussed  in  the  main 
text,  analytical  results  for  these  moments  cannot  be  derived  when  the  full 
Monin-Obukhov  form  for  the  dimensionless  wind  gradient  (eq  (15))  is  used. 
Therefore,  in  this  appendix,  I  consider  two  limiting  approximations  for  the 
local  wind  and  temperature  profiles,  based  on  the  limiting  cases  of  neutral 
and  free-convection  stability. 

In  the  neutral  limit,  'I'm  (C)  =0  the  local  velocity  has  a  logarithmic 
profile: 


U  T  Z 

Ul  (z)  =  [{U  (z))  +  (z)]  +  eyGy  (z)  =  (e^;  cos  0  +  ey  sin  6)  —  In  —  .  (B-1) 

K  Zq 

Evaluating  this  equation  at  z  =  Zr  and  taking  the  magnitude,  one  has 

u,L  =  Kln-^-  J {Ur  +  G,,f  +  Gl  .  (B-2) 

Now,  substituting  equations  (B-2),  (29),  and  (8)  into  (9)  results  in 


cos™  0  sin"  0)  = 


/oo  poo  p  -| 

/  {Ur  -|-  Gx)'^  -|-  Gy 
-oo  J  —OO  ^ 


{q—m—n)/2 


In^  {Zr/Zo)  I 

( 

x([/.  +  U.)™G^xp 


dGx  dGy  . 


Converting  to  cylindrical  coordinates  p  and  6  such  that  pcosO  = 
{Gx  +  Ur)  / and  psinO  =  Gy//3w^:,  the  double  integral  becomes 

—2  r°°  /■27r  _ 

,  cos™  0  sin”  9)  =  — - — -e~^^  /  /  cos™  6  sin”  9  exp  {—p^  -|-  2Ur  cos  9)  d9  dp  , 

‘  Trln^  (2;,./zo)  Jo  Jo 

(B-4) 

where  u^l  =  u^l/Pui*  and  Ur  =  Ur/fdw^. 

We  can  perform  the  integral  over  9  using  the  integral  representation  of  the 
modified  Bessel  function  (e.g.,  Abramowitz  and  Stegun  (1965),  eq  (9.6.19)): 

In{z)  =  —  [  COS  {n9)  d9  .  (B-5) 

^  Jo 
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Considering  specific  values  of  m  and  n,  we  therefore  have 

{<l) 

{uIl  cos  e) 

{uIl  COS^  e) 

{ulisiT^e) 

The  latter  two  relationships  follow  from  the  identities  2  cos^  0  =  1+  cos  26 
and  2  sin^  6  =  1—  cos  26. 

The  remaining  integrals  in  equations  (B-6)  through  (B-9)  all  have  the  form 

I j  (bt)  dt  . 

The  strategy  for  performing  these  integrals  depends  on  whether  k  —  j  is 
even  or  odd.  When  it  is  odd,  the  starting  point  is  the  following  relationship, 
which  applies  to  k  —  j  =  1: 

rOD  VI/ 

m  dt  =  .  (B-10) 

This  relationship  can  be  found  in  standard  tables  of  integrals.  For  exam¬ 
ple,  it  is  equivalent  to  equation  (11.4.29)  in  Abramowitz  and  Stegun  (1965), 
after  applying  their  equation  (9.6.3).  The  cases  k  —  j  =  3, 5, . . .  can  be  re¬ 
duced  to  equation  (B-10)  by  repeated  applications  of  the  recurrence  formula 
In+i  (x)  =  In-1  (x)  —  (2n/x)  In  (x).  Solutions  for  even  k  —  j  are  based  on  the 
following  result  (Abramowitz  and  Stegun's  eq  (11.4.32))  for  k  —  j  =  0: 

m  dt  =  (I)  •  (B-11) 

The  cases  k  —  j  =  2,4,...  can  be  determined  through  integration  by  parts. 
In  the  formula  f  udv  =  uv  —  f  v  du,  repeatedly  set  v  =  until  integrals 
having  the  same  form  as  equation  (B-11)  result.  For  reference,  the  special 
cases  of  interest  are  listed  here: 


2k-'?  —2  r°°  9  _ 

=  '  I  '‘P  ■  (B-6) 

_ 2  9  _  _ 

=  V+l — [loijlUrp)  +  l2{j2.Urp)]  dp  ,  and  (B-8) 
In  {Zr!  Zq)  Jq 

_ 2  9  _  _ 

=  llPH+I+f"'  i  (2U.P)]  dp  .  (B-9) 
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{u^l) 
{u^L  COS  9) 
{uIl  cose) 
{uIl  COS^e) 

{uIl  sin^  e) 


{bt)  dt 
“*^/o  (bt)  dt 
(bt)  dt 
“*^/o  (bt)  dt 
~^^Ii  (bt)  dt 
h  (bt)  dt 


(B-12) 

(B-13) 

(B-14) 

(B-15) 

(B-16) 

(B-17) 


With  these  integrals,  we  can  now  write  down  the  following  equations  for 
moments  of  the  local  friction  velocity: 


\/^ln(^r/^o)„-f7!/2r 
- e  Jo 

K 


2\V?  {Zr/zo) 


(B-18) 
(B-19) 
(B-20) 
,  (B-21) 
(B-22) 
(B-23) 


Since  Iq  (0)  =  1  and  Ii  (0)  =  0,  in  the  limit  Ur  0,  these  equations  become 

^ln(2;,./2:o) 


{kD  = 

{u*l)  = 


K 

TTK 


21n(2;r/2;o) 

{u^L  cos  9)  =  {uh  cos  0)  =  0  ,  and 


( 


uIl  cos' 


u^L  Sin 


'9)  = 


21n^  {zr/zo) 


(B-24) 

(B-25) 

(B-26) 

(B-27) 


The  analysis  for  the  free-convection  limit  proceeds  similarly  to  the  neu¬ 
tral  limit.  Instead  of  equation  (B-1),  we  use  (C)  =  InC  +  3  (— 
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(Grachev  et  al,  1998),  resulting  in 


\J L  {z)  —  Gj;  [{U  (z))  +  Gx  (2:)]  +  GyGy  (z) 


=  (e^;  cos  6  +  Gy  sin  9) 


-I 

ZoJ 


\  1/3 


Zi 


1/3' 


(B-28) 


instead  of  equation  (B-1).  Evaluating  this  equation  at  z  =  Zr  and  taking  the 
magnitude,  we  find  that 


(B-29) 

Therefore 


.4/3„V3\<?/2 


COO  /*27r 


..  \  1/3  / X 1/31  -5/2 


iV' 

oj  \ZrJ 


roo  pZn 

/  /  /9'^^^’'"^cos™'0sin’^0exp  (— /)^  +  2[/rCos0)  . 

Jo  Jo 


(B-30) 

(B-31) 


The  integral  over  6  proceeds  in  the  same  manner  as  the  neutral  limit.  When 
q  is  even,  we  can  also  determine  the  remaining  integral  over  p  using  the 
previous  methodology.  Integration  for  odd  q  leads  to  confluent  hypergeo¬ 
metric  functions.  However,  since  the  free-convection  limit  should  be  useful 
only  for  small  Ur,  we  may  simplify  the  integrations  by  setting  Ur  =  0.  With 
this  approximation. 


1  /  4/3  1/3  \  1/2 

(s’^cos^Osin”®)  =  -  j 


ZoJ  V^r 


1/3' 


-q/2 


roo  r27T 


'  rZTT 

/  cos”^  6  sin”  6  exp  (— p^)  d9  dp  . 

/o  Jo 

The  integrations  now  yield 


{<l) 

{ul^cos‘^9) 


{Kl  cos 


=  r((?/4  +  l) 

V 

=  0, 

=  (ii*Lsiii^0>  = 


4/3„y3\  9/2 


K  '"a 


M 


3/3 


ZoJ 


1/3 


1/31  -9/2 


'(d/4  +  1)  ( 


3/3 


zo. 


)"-(!) 


1/3' 


1  -9/2 


(B-32) 

(B-33) 


(B-34) 

(B-35) 

(B-36) 
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and 


{ul^sin^e) 


r  {q/4  +  1)  / 

2  3/3 


(B-37) 


Figures  B-1  and  B-2  compare  the  neutral  and  free-convection  limits  to 
full  Monin-Obukhov  similarity  for  moments  involving  arvd  cos  6.  The 
plots  show  that  neither  the  neutral  nor  the  free-convection  limits  provide 
a  very  satisfactory  approximation.  As  Ur  increases,  the  neutral  limit  pre¬ 
diction  converges  very  slowly  to  full  Monin-Obukhov  similarity,  differing 
by  more  than  10  percent  even  at  Ur  =  10.  The  free-convection  prediction 
very  roughly  approximates  (to  within  about  30%)  the  moments  and 
(ul]^  cos^  6)  when  [7^  <C  1. 


Figure  B-1.  Statistical 
moments  (u*l)  and 
(m*l  cos  9)  as  a  function 
of  normalized  wind 
speed  Ur  -  Local 
equilibrium  predictions 
that  use  full 
Monin-Obukhov 
similarity,  as  well  as 
approximate 
predictions  in  neutral 
and  free-convection 
limits,  are  shown. 
Roughness  length  is 

(Free-convection 
prediction  for 
{ut-L  COS  9)  is  zero.) 


Figure  B-2.  Same  as 
figure  B-1,  except  for 
statistical  moments 

(u*L  cos^  0)  and 

(uIl  cos^  ^)- 

(Free-convection 
prediction  for 
(uli^  cos  0)  is  zero.) 
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